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ABSTRACT 

We study the estimators of various second-order weak lensing statistics such as the 
shear correlation functions £± and the aperture mass dispersion (M£ p ) which can 
directly be constructed from weak lensing shear maps. We compare the efficiency with 
which these estimators can be used to constrain cosmological parameters. To this end 
we introduce the Karhunen-Loeve (KL) eigenmode analysis techniques for weak lensing 
surveys. These tools are shown to be very effective as a diagnostics for optimising 
survey strategies. The usefulness of these tools to study the effect of angular binning, 
the depth and width of the survey and noise contributions due to intrinsic ellipticities 
and number density of source galaxies on the estimation of cosmological parameters 
is demonstrated. Results from independent analysis of various parameters and joint 
estimations are compared. We also study how degeneracies among various cosmological 
and survey parameters affect the eigenmodes associated with these parameters. 

Key words: Cosmology: theory - gravitational lensing - large-scale structure of 
Universe - Methods: analytical, statistical, numerical 



1 INTRODUCTION 

Weak gravitational lensing studies (for a review see e.g. Mel- 
lier 1999; Bartelmann & Schneider 2001; Refregier 2003; van 
Waerbeke & Mellier 2003) probe the correlation of the ob- 
served ellipticities of background galaxies. Cosmic shear pro- 
vides unbiased estimates of the matter power spectrum and 
hence can help us to place independent constraints on cos- 
mological parameters (see e.g. Hu 1999) . The parameter de- 
generacies from weak lensing surveys are very different from 
those present in CMB studies (van Waerbeke et al. 2002). 
This means that even most precise measurements from CMB 
observations can be improved by joint analyses with weak 
lensing data (Contaldi et al. 2003; Hu & Tegmark 1999). 
Thus, the recent detection of weak lensing by the large- 
scale structure (Bacon et al. 2000; Kaiser, Wilson & Luppino 
2000; Wittman et al. 2000.; van Waerbeke et al. 2000, 2002; 
Maoli et al. 2001; Refregier, Rhodes & Growth 2002) has 
opened a completely new window to the Universe. While 
first generations of cosmic shear surveys have demonstrated 
the feasibility of weak lensing studies in constraining the 
dark matter power spectrum parametrised by Q m , as and 
the shape parameter F (Refregier 2003; van Waerbeke, Mel- 
lier & Hoekstra 2005), future surveys will be able to probe 
much larger scales and therefore allow us to study the linear 
regime directly and put more stringent bounds on cosmolog- 



ical parameters such as the equation of state of dark energy 
and its time variation. 

As survey sizes grow the amount of data that is to 
be processed also increases. It is then essential to study 
data compression techniques which can effectively compress 
the data vector to a lower-dimensional space. The effect of 
such a compression on the Fisher information matrix can be 
studied by analysing the Karhunen-Loeve (KL) eigenmodes 
(Karhunen 1947; Loeve 1948) associated with the Fisher 
matrix. These techniques are already in use in cosmic mi- 
crowave background studies (see e.g. Bond 1995 for COBE 
& FIRAS data, Bunn 1995 for COBE data and Tegmark, 
Taylor & Heavens 1997, hereafter TTH, for a general review 
of the KL eigenmode analysis). In the context of galaxy 
surveys, Matsubara, Szalay & Landy (2000) applied KL 
methods to the Las Campanas redshift survey, Hamilton, 
Tegmark & Padmanabhan (2000) studied PSCz using KL 
techniques and Szalay et al. (2003) have used this method 
to study data from SDSS. Watkins et al. (2001) have used 
these techniques to analyse the peculiar velocity surveys for 
removal of small-scale, non-linear modes from large-scale 
linear modes which retain cosmological information. Other 
studies using the KL eigenmode methods for velocity field 
surveys include Hoffmann & Zaroubi (2000) and Silberman 
et al. (2001). Recent studies by Huterer & White (2005) 
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have underlined its importance for filtering unwanted non- 
cosmological information at small angular scales. 

The data vector in weak lensing surveys which is used 
to extract cosmological information is usually the shear cor- 
relation at various angular separations, estimated from the 
shapes of the observed background galaxies. Thus, it rep- 
resents an already compressed data set for which likeli- 
hood analyses are feasible although further compression can 
speed up the analysis. Moreover, for third-order statistics, 
the number of measured triangle configurations can be very 
large even in a binned way and and further compression of 
the data might be necessary to place constraints in a high- 
dimensional cosmological parameter space. Therefore, the 
study of KL methods in the context of cosmic shear will be 
very useful for future weak lensing surveys. 

Being not only an effective tool for data compression, 
the KL eigenmode analysis can also help to understand the 
specific linear combination of angular scales and the red- 
shift range that contributes most significantly to a specific 
estimator. A KL eigenmode analysis of weak lensing observ- 
ables in redshift space has already been done in by Heavens 
(2003), where the distribution of observed ellipticities was 
directly used as data vector. We extend these analyses to 
second-order shear statistics as functions of projected an- 
gular scales and for realistic survey geometries. The latter 
issue is particularly interesting as sky coverage in any weak 
lensing survey will always be non-trivial due to the presence 
of bright foreground stars and galaxies. Our study therefore 
gives us useful clues about the survey design. We provide a 
KL eigenmode analysis for the second-order shear measures 
?±j ftot and (Mf„) (to be defined later). The covariance 
of these estimators were studied in detail by Schneider et 
al. (2002) to assess the effect of source ellipticity dispersion 
and finite sky coverage on the estimation of cosmological 
parameters. The influence of survey geometry was investi- 
gated in Kilbinger & Schneider (2004). We use those results 
to construct the Fisher matrix to forecast parameter con- 
straints. In contrast to those previous studies, where the 
above mentioned second-order shear statistics were used to 
determine cosmological parameters, we also consider the pa- 
rameter dependence of the covariance of these estimators 
and its influence on the Fisher matrix. Moreover, the KL 
approach allows us to study in detail how the convergence 
power spectrum is sampled on different scales by a given 
shear survey. 

The paper is arranged in the following way. Section Sj5| 
is devoted to a recapitulation of the basics of the KL eigen- 
mode analysis, reformatted in the weak lensing context. In 
section §3, we describe our simulations and survey strategies 
in detail (§3.1 - §3.3) and present the results for the eigen- 
mode analysis (§3.4), including the effect of survey strategy, 
noise and binning. We discuss our findings in §4. 



2 NOTATIONS AND FORMALISM 
2.1 Survey characteristics 

We consider several different survey strategies in this work. 
First, individual images are distributed randomly (but non- 
overlapping) within circular patches on the sky, where 'im- 
age' means one individual field of view of size 13' x 13'. The 



individual patches are uncorrelated, thus the largest scale 
where the shear correlation can be measured is twice the 
patch radius. The random image placement accounts for the 
fact that for realistic surveys, bright stars and foreground 
galaxies should be avoided and therefore sparse sampling of 
a sky region is necessary. 

A survey consists of P patches of radius R, each patch 
containing N images. The total number of images is 300 for 
all patch geometries, corresponding to 14.1 square degrees 
for the entire survey. 

We use N = 30 and R = 100', 140', respectively. Note 
that the number of patches in a survey is P — 300/iV. We 
denote these patch geometries by the two numbers (N,R), 
thus we have the configurations (30, 100') and (30, 140'). 

Second, we define two surveys which consist of single, 
uncorrelated fields of view with the same total area of 14.1 
square degrees. The first survey consists of 75 uncorrelated 
26' x 26' square images - typical fields of view for current 
wide- field cameras like WFI on the ESO 2.2-m telescope. 
The second survey represents 12 independent 65' x 65'- 
fields, corresponding to new-generation wide-field cameras 
like MegaCam/CFHT. We denote these surveys with 75-26' 2 
and 12 • 65' 2 , respectively. A sketch of the geometry of the 
surveys used in this work can be found in Fig. 

If not indicated otherwise, the number density of source 
galaxies is n ga i = 30 arcmin -2 . This number density of high- 
redshift galaxies which are usable for weak lensing shape 
measurements can be achieved with high-quality ground- 
based imaging data from a 4 m-class telescope. The source 
galaxy ellipticity dispersion is o> = 0.3, if not stated other- 
wise. 



2.2 Second-order shear statistics 

The light of distant galaxies is deflected while propagat- 
ing through the tidal gravitational potential of interven- 
ing matter inhomogeneities which constitute the large-scale 
structure in the Universe. The image of a galaxy gets dis- 
torted, and in the weak lensing approximation, the relation 
e — e 8 +7 holds between the observed ellipticity £ of a galaxy, 
its intrinsic ellipticity e s and the shear 7 = 71-H72. From the 
shear, one can in principle reconstruct the projected matter 
density or convergence «. 

For two points separated by a vector 6 with polar angle 
tp, one defines the tangential and cross-component of the 
shear as 

7t = - 5R ^e~ 2lv ^ and 7 X = - 9; (7e~ 21 ^) • (1) 

The two shear two-point correlation functions (2PCF) and 
their relation to the convergence power spectrum P K are 
(Kaiser 1992) 

OO 

£±(0) = (7t7t) ± <7x7x> = ^ J d££P K (£)J 0A (m, (2) 



where the first-kind Bessel function Jo (J4) corresponds to 
the '+' ('— ') correlation function. 

A very useful second-order statistics is the dispersion of 
the so-called aperture mass (Kaiser 1995, Schneider 1996), 
which is the weighted tangential shear integrated over a cir- 
cular aperture. When using a polynomial compensated filter 
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Figure 1. Sketch of the two kinds of survey strategies used in 
this work. On the left, a patch of 100/ radius is filled sparsely with 
30 images of size 13', a survey consists of 10 such patches. For 
comparison, on the right the field sizes of surveys consisting of 75 
uncorrelated 26' X 26' images and 12 independent 65' X 65' fields, 
respectively, are shown. 



function as weight (Schneider 1996, Schneider et al. 1998), 
the dispersion of the aperture mass is related to the power 
spectrum by 



oc 

«(*)> = £ / 



df.£P K (£) 



24 J 4 (£0) 

(ley 2 



(3) 



The aperture mass statistics is a measure of the E-mode 
only of cosmic shear. A similar statistics {Mj_) can be de- 
fined as the weighted cross-component of shear in an aper- 
ture, which is a measure of the B-mode only. The aperture 
mass statistics yields the separate measurement of the E- 
and B-mode which allows one to quantify contaminations 
to the gravitational shear signal. A B-mode mainly comes 
from systematic measurement errors like imperfect PSF cor- 
rection and/or intrinsic alignment of source galaxies, whose 
shape correlations are to be measured. 

The 2PCF are easily obtained from real data however 
complicated the survey geometry might be. The two aper- 
ture statistics can be expressed in terms of the 2PCF: 



Add 
~29 T 



,(4) 



where T+ and T_ are given explicitly in Schneider, van 
Waerbeke & Mellier (2002). In the absence of a B-mode, 
both terms on the right-hand side of the previous equation 
are of the same amplitude, resulting in (Mj_) = 0. 

An unbiased estimator of the shear 2PCF J5J for a set 
of (logarithmic) angular bins is given by 

1 



e±(*) = 



(5) 



where the sum goes over all pairs of galaxies, whose angular 
separation is in the angular bin corresponding to e< is the 
(complex) ellipticity of a galaxy at position 9% and N v {$) is 
the number of pairs in the #-bin. Galaxies measured from 
real data might carry weight factors reflecting the individual 
noise and precision in the shape measurement. Using these 
weights, one can define a more optimal estimator for the 



2PCF; however, we do not include this effect and assume 
equal weights for each galaxy. 

Unbiased estimators for the dispersions of the two aper- 
ture mass statistics (Mf p y and (M±) respectively, are 



)T- 



(t) 



(6) 



where the limit I of the sum must be chosen such that 9 is 
half the upper limit of the 7-th bin. For simplicity of nota- 
tion, we will from now on identify the estimators with their 
corresponding shear statistics, and write and (Aff p ) 

for the estimators. 

These second-order shear estimators are used exten- 
sively in the literature, including independent studies in- 
volving £± or the joint data vector £ t ot = (£+,£-)• We do 
not take into account the shear dispersion in an aperture 
(|7| 2 ). Although this estimator has a high signal-to-noise 
ratio and has been used extensively in the literature, it does 
not separate the E- from the B-mode, unlike the aperture 
mass statistics. 

Schneider et al. (2002) calculated the covariance matri- 
ces C of the estimators defined in this section, in the case 
of a Gaussian shear field. The covariances of the 2PCF con- 
sist of three terms: A pure shot noise term (D), originating 
from the dispersion of the intrinsic ellipticities of the source 
galaxies, present only on the diagonal, a cosmic variance 
term (V), and a mixed term (M): 



C++ = (#+,01; = D + M H 

C__ =C(£_,0i;|_,aO = D + M- 
C+- =C(£h*i; = M H 



+ + v, 

+ V 

+ 



(7) 



The individual terms in full length are given in Schneider 
et al. (2002). The covariance of both of the aperture mass 
statistics, C(A4±), can be calculated from the 2PCF covari- 
ance matrices using eq. ©. 

Derivatives of the covariances with respect to cosmo- 
logical parameter are easily obtained by differentiating eqs. 
Q. The shot- noise term, being independent of cosmology, 
does not contribute, and the cosmic variance term, being 
quadratic in the correlation function, has a larger influence 
than the mixed term. 

The covariances of these estimators were studied in de- 
tail by Schneider et al. (2002), where covariance matrices 
were constructed using an ensemble average over galaxy 
positions. Kilbinger & Schneider (2004) calculated the co- 
variances directly by summing over a specific realisation of 
galaxy positions which we follow here. This has the advan- 
tage of being able to handle the discreteness resulting from 
the finite galaxy number density in a correct way. More- 
over, incorporations of realistic survey strategies is done in 
a natural way. 

We note here that the cosmic variance term V only con- 
tains the Gaussian or unconnected part of the four-point 
correlator of shear. On scales below ~ 10 arc minutes, the 
non-Gaussianity of the shear field gets important (van Waer- 
beke et al. 2002; Scoccimarro, Zaldarriaga & Hui 1999). On 
scales below ~ 1 arc minute the shot noise term D domi- 
nates over V (Kilbinger & Schneider 2004), thus under the 
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Gaussian assumption we expect to slightly underestimate 
the covariances in this transition regime between 1 and 10 
arc minutes. 

2.3 Fisher information matrix 

The Fisher matrix (Kendall & Stuart 1969; TTH) is denned 

as 

Fal3 \ dp a dpp } ( dp a dpp ) p= p ^ 

where the likelihood C depends on a vector of model pa- 
rameters p = (pi,p2, •■•) and p denotes the true parameter 
values. The Fisher matrix is the expectation value of the 
Hessian matrix of (— ln£) at p — p and is roughly spoken 
a measure of how fast C falls off from its maximum. The 
smallest possible variance Ap a of any unbiased estimator of 
some parameter p a is given by the Cramer-Rao inequality 

Ap a > y/{F-^) aa - (9) 

the right-hand side of this inequality is called the minimum 
variance bound (MVB). If only one parameter is to be de- 
termined from the data and all others are fixed, the MVB 
simplifies to (F aa ) _1 ' 2 . 

In the case of a Gaussian probability distribution func- 
tion for the parameters, the Fisher matrix becomes (Bunn 
1995; Vogley & Szalay 1996; TTH) 

F a0 = itrfC^C.c C _1 C^ + C^M^] . (10) 

C, a = dC/dpc denotes the derivative of the covariance with 
respect to the a th parameter. M a p = p a p t /3 +fJ, /3 p t a where 
p, is mean of the data vector x, (x) — p. In our case, the 
n- dimensional data vector x consist of one or a combination 
of the second-order statistics of shear as defined in Sect. 12.21 
(see also Sect. I3.l1 for more details). 

The m = 7 parameters p a (cosmological parameters 
plus source galaxy redshift distribution characteristics) 
which we consider here are given in Sect. 13.21 

2.4 Karhunen-Lofeve (KL) eigenmodes and 
eigenvalues 

A general linear data compression can be written as 

x = Tx, (11) 

where from the n-dimensional data vector x a new h- 
dimensional data set x is constructed via the h x n ma- 
trix T. This means that the expectation value (x) of the 
original data vector x transforms to (x) — T{x) and hence 
the covariance C = (xx l ) — {x)(x) t changes accordingly to 
C = TCT'. In our case, X consists of the shear correlation 
functions and £_ or the aperture mass dispersion (M^ p ), 
representing n data points measured for various angular 
scales, see Sect. 13.11 An example for a data compression 
is already given in eq. @: (Aff p ) is a linearly compressed 
version of £+ and £_ , the compression matrix T depends in 
this case on the functions T+ and T_ . 

The Fisher matrix F corresponding to the original data 
vector x changes to F, given by 

F a/3 = itr[(TCr)- 1 (TC,„T t )(TCT t )- 1 (TC, /J r) 



+(TCT t )- 1 (TM Q/3 T t )] (12) 

Clearly, when the dimensionality n and h of the original and 
the transformed data vectors are the same, the transforma- 
tion is a similarity transformation and the Fisher matrix 
remains unchanged, F = F. 

The goal of the KL analysis is to find a compression 
matrix T for h < n without loosing too much of information 
and to estimate cosmological parameters with as small error 
bars as possible. Considering the minimum variance bound 
of some cosmological parameter p a , we seek to maximise 
the Fisher matrix element F aa . 

Following TTH, we first consider the two simple cases of 
a constant mean fi and a constant covariance C, where con- 
stant means with respect to the parameter vector p. Then, 
we optimise for several parameters simultaneously with the 
aid of a singular value decomposition (SVD). 

For cosmic shear, the assumption of a constant covari- 
ance seems quite natural, since in most cases of weak lensing 
observations, the covariance is not calculated analytically 
but extracted directly from the data (e.g. van Waerbeke et 
al. 2002) and thus it is constant and not depending on cos- 
mological parameters. However, the situation is different for 
the mean, since with a model parameter independent mean, 
it is not possible to even define a likelihood function. And al- 
though within the Fisher matrix formalism, a constant mean 
can be considered in a consistent way, this mean first has to 
be found using e.g. the maximum likelihood estimator. But 
the distinction into a constant mean and a constant covari- 
ance case is of rather technical nature, and the general case 
of parameter dependend mean and covariance can be consid- 
ered by combining the two, see Sect. 12. 431 Thus, a strategy 
to obtain cosmological parameters and their minimum vari- 
ance bound from data would be first to define a likelihood 
function C using a parameter dependent mean (and param- 
eter dependent covariance if desired) to find the maximum 
likelihood parameter p . Second, the Fisher matrix is calcu- 
lated by evaluating the second derivatives of C at the point 
p = p . A KL analysis and data compression can be under- 
taken by first considering the two cases of constant mean and 
covariance independently, and later combining both cases. 

Note that in our case, the "mean" is a second-order 
statistics (the shear correlation), and the covariance is of 
fourth order (although reduced to only depend on second- 
order because of the assumption of a Gaussian shear field). 
We therefore choose a different approach than e.g. Heavens 
(2003) or Seljak (1998) who directly use the galaxy elliptic- 
ities as data vector. In their case, the (zero) mean is of first 
order, and the covariance contains the information about 
the power spectrum. 

2.4-1 Constant Mean 

In the case where the mean p, is known and constant, the 
second term in 112t vanishes, since p, a = and therefore 
M Q /3 vanishes, too. We first consider a compression matrix 
which consists of only one row vector, T = b . The task 
of maximising the Fisher matrix diagonal element F aa in 
order to minimise the error bar on the cosmological param- 
eter p a is then equivalent to maximising the eigenvalue A 
of the generalised eigenproblem C, Q 6 = ACfo. By Cholesky- 
decomposing the symmetric and positive definite covariance 
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matrix C = LL*, it can be reduced to the ordinary eigen- 
problem 

(L t C, a L- t )L t b = xVb. (13) 

Solving this equation for all n orthogonal eigenvectors L'b^ 
gives us n real eigenvalues Xk- A compression matrix T is 
then constructed containing as row vectors the first h eigen- 
vectors which have been sorted by the absolute value of their 
corresponding eigenvalues. The matrix T represents a set of 
eigenvectors, rank-ordered according to their signal-to-noise 
ratio. 

The KL eigenmodes defined in this way satisfy an or- 
thogonality relation. Eigenmodes with high eigenvalues or 
low rank numbers contain more information about a specific 
parameter, whereas the ones with small (absolute) eigenval- 
ues and high rank numbers contain almost no additional in- 
formation. As will be shown in Sect. 13.31 we typically achieve 
a compression by a factor of nearly two for the independent 
analysis of various parameters with getting the same MVB 
as for the uncompressed case. 

An individual eigenmode contributes to the mea- 
surement error for the parameter p a as 5p a — l/|Afc|. 
Thus, the signal-to-noise ratio for this mode therefore is 

PaiSpa)' 1 = Pa|Afc|. 

The new Fisher matrix diagonal element F aa corre- 
sponding to a KL compression with h eigenmodes is simply 
(TTH) 

n 

F aa = ^X 2 k . (14) 
fc=i 

In our analysis (Sect. l3~3"l . we repeat the data compression, 
starting from the second-order statistics £± and (M% p ) given 
for n angular separations and plot for various cosmologi- 
cal parameters p a the associated error Ap a = lf(Faa) 1 ^ 2 
as a function of h, where h < n is the dimension of the 
compressed data vector x, see Fig. [5] Although being a de- 
creasing function of h, the error reaches a constant plateau 
for some no < n, thus the original error for the uncom- 
pressed case is recovered before all KL modes are used for 
the parameter estimation. We find in most of the cases that 
a plateau is reached for ho n/2, thus a compression fac- 
tor of nearly two is possible without any loss of information. 
This compression is with respect to the number of angular 
separations n at which the second-order shear statistics are 
measured, representing a binned version of the data vector 
comprising of all observed galaxy pairs. Note that the com- 
pression factor depends of course on the original number of 
bins — we comment on the binning in Sect. I3.3"H1 

Since in the case of a constant mean, the Fisher matrix 
contains products of the inverse covariance and the (deriva- 
tive of the) covariance, it is independent of the survey area, 
and very little sensitive of the survey geometry. 



2.4-2 Constant covariance 

In contrast to the constant mean case, there is only one 
eigenvector which contains all of the available information 
when the covariance is known and independent of the pa- 
rameters. In this case, the first term in 1121 vanishes, and 
the corresponding eigenproblem is 



(L~ 1 M aQ L _t )L t 6 = AL'fo (15) 

which has only one non-trivial solution l^bo = L 1 fi a with 
eigenvalue Ao = 2 1 L — 1 /x | = tr [C _1 M Qa ] . 

Inserting the only eigenvector bo = C _1 /x a for the 
compression matrix T into 11 'it . we find that the modified 
Fisher matrix is the same as the original one, F a p = F a p = 
/x t a C _1 /x g. This is not surprising since all data is collected 
into one mode, and no information is lost. 

If geometrical effects of the survey are neglected, the 
covariance is proportional to one over the survey area A. 
Thus, in the case of constant covariance the Fisher matrix 
is proportional to A, in contrast to the constant mean case, 
where F is independent of A. For reasonable large cosmic 
shear surveys, the second term in HOH is therefore dominant 
over the first one fSect. |3~3.5H . 

2-4-3 General case 

TTH describe how the general case (when neither mean nor 
covariance is constant) can be treated efficiently, by simply 
adding the one eigenmode from the constant covariance case 
to the h modes from the constant mean analysis. However, as 
mentioned in the previous section, the constant covariance 
eigenmode is (for reasonably large survey areas) dominant 
over all the other modes and contains the bulk part of the 
information about cosmological parameters. Thus, we do not 
consider to combine these two cases. 



2-4-4 Joint parameter estimation 

For the independent estimation of a parameter p a , the KL 
method is devised to minimise the associated error bar 
Ap a = (Faa) -1 / 2 in a suitably rotated basis. However, for 
the case of joint parameter estimation, the object to be min- 
imised is the diagonal of (F 1 ) 1 / 2 which is a more demanding 
optimisation problem. We follow TTH who discuss an alter- 
native approximate technique which virtually does the same. 
The individual compression matrices T a , optimised for the 
independent estimation of the parameter p a , are arranged 
into a new matrix after multiplying each row (which corre- 
sponds to an eigenvector) with its corresponding eigenvalue 
as T = (AjTi, . . . ,A m T m ). Here, A a = diag(A a i, . . . ,X a n) 
is a diagonal matrix containing the eigenvalues correspond- 
ing to the individual eigenanalysis of the a th parameter. 
T will contain a lot of redundant information, because of 
unavoidable near-degeneracies between parameters, e.g. be- 
tween fl m and as or between F and n s . 

In order to separate useful from redundant information, 
this new matrix T is factorised using a singular value decom- 
position (SVD), 

nm n n nm 

n{ ( Axil | •■■ | A ro T* m ) = ( U ) ■ ( £ ) • ( V 1 ) }„ 

where 11*1) = V*V = 1 and S = diag(<7i) is a diagonal ma- 
trix containing the singular values <7;, which can be inter- 
preted as generalised eigenvalues. After sorting the columns 
of U* according to the absolute values of the corresponding 
a the final compression matrix is Tj i nt = U'. Modes with 
high singular values contain the bulk information about the 
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Figure 2. Ap a = (Fan)" 1 ' 2 as a function of the mode num- 
ber h is plotted, according to 1141 for the constant mean case 
and various cosmological parameters. Solid lines correspond to 
the independent parameter estimation (Sect. ~T.4.H . dashed lines 
represent the joint parameter analysis (Sect. ~~4.4l , The four pan- 
els correspond to the shear statistics f+, their combination 
ftot and (Aff ) (see Sect. l~~Tl The survey geometry is (30, 100'). 
The maximum mode number n equals the number of bins, which 
is 38 for ftot, 19 for both £+ and §_ and 14 for (Mf p ). 

cosmological parameters, whereas row vectors of Tj i nt cor- 
responding to vanishing or very small m capture redundant 
or almost redundant information. Inserting Tjoint into (1121) . 
the Fisher matrix after the final compression is easily calcu- 
lated. 

If n a is the number of modes that were used in the 
optimisation process for the a th parameter, the joint matrix 
Tj i n t has fi — h a columns. Based on the amplitude 

of the singular values, one can choose a final mode number 
Wjoint < n' fixing the compression factor. As an example, 
the transformed diagonal elements of the Fisher matrix are 
plotted as a function of fijoint in Fig. [5] where the saturation 
plateau can be used to determine the final mode number. 

Also for the constant covariance case (Sect l2.4"~l . a joint 
parameter estimation analysis can be done. We combine the 
m eigenvectors resulting from the m individual parameter 
analyses to the nxm-matrix T, after having multiplied them 
by their respective eigenvalues. The singular value decom- 
position of T will result in m singular values, which give us 
information of the degree of degeneracy between the param- 
eters. We will discuss results for this analysis in Sect. 13.31 

2.4-5 Window functions 

The linear data compression defined in 11 U generates a new 
data vector x which in the case of single parameter optimisa- 
tion represents a set of uncorrelated unit-variance Gaussian 
random variables (TTH). The components Xi of the new 



data vector contain pairwise independent and uncorrelated 
information about cosmology. In our case the data consists 
of second-order shear statistics and the cosmological infor- 
mation is contained in the convergence power spectrum. 

Using the KL eigenmode technique, we can study in 
detail how the power spectrum is sampled in order to yield 
uncorrelated data points. By comparing modes correspond- 
ing to high (low) eigenvalues, we can see which scales carry 
much (little) of the cosmological information. 

For any second-order statistics, the dependence on the 
power spectrum is encoded in a window function. For the 
two shear 2PCF £+ and , these are the broad-band Bessel 
functions Jo and J4, respectively, see eq. J~j; the filter func- 
tion for the aperture mass statistics (Mf p ) is the narrow- 
peaked function [24 J 4 (»?)/?? 2 ] 2 , see ~~). 

For each component Xi of the new data vector, we de- 
fine a new window function Wi, which is a combination of 
the original filter functions associated to different smoothing 
scales 8j. In the case of the 2PCF, the new data vector is 

1 f°° 

If = T^ ± {6 3 ) = —J d££P K (£)wt(£); 

W?{£) = 5^Z£ J °.*(#J)- (16) 

3 

Similarly, for (M 2 p ) we obtain the new window functions 
Wf{£) = Y,- Tij[24 J 4 (^)/(^) 2 ] 2 . In the case of the com- 
bined vector f tot = the window functions Wl° % are 
linear combinations of and W~ . 

In case of a constant mean, the rows of the matrix T 
contain the n transposed eigenvectors b, see Sect. l2.4~T1 For 
the constant covariance case, there is only one eigenmode, in 
this case To = /x t a C _t . We denote all quantities correspond- 
ing to the constant covariance case with a subscript '0', e.g. 
£0 > Wo" e ^ c - F° r t ne joint parameter estimation (Sect. l2~4.4ll . 
window functions are defined analogously, by inserting the 
joint matrix Tjoint into cq. 1161) . 

The study of these window functions, which is presented 
in the next section, will provide us with insights about how 
the convergence power spectrum is sampled in order to con- 
strain cosmological parameters with cosmic shear. It allows 
us to investigate how different Fourier modes are probed by 
a given survey strategy. 



3 NUMERICAL RESULTS 

The impact of various survey strategies, noise structures and 
the effect of binning are considered in the determination of 
the eigenmodes. We study separately the constant mean case 
(see Sect. I2.4~TI and the constant covariance case (see Sect. 
12.4.2^ . Eigenmodes associated with the independent analy- 
sis of a single parameter and the joint analysis of several 
parameters (see Sect. I2.4~~t are investigated. 

3.1 The input data 

Our input data vector x is one of one the following second- 
order statistics of cosmic shear, as defined in Sect. 12.21 

1. The correlation function measured at various 

angular scales The covariance is denoted as C++. 
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2. The correlation function measured at various 
angular scales Qi. The covariance is denoted as C 

3. The joint combination ftot of both data vectors f + and 
£_ where the joint covariance matrix Ctot is constructed 

out of blocks C++, C on the diagonal and C+_ and its 

transpose on the off-diagonal. 

4. The aperture mass statistics (Mf p (#j)) for various 
aperture radii 9i and its covariance C(M + ). 

These statistics are calculated from the convergence 
power spectrum using eqs. J2J and 0. For the 2PCF, we 
use 19 angular logarithmic bins, the smallest separation be- 
tween two galaxies considered being 0.5 arcmin. Since the 
calculation of the covariances and their derivations are very 
time-consuming, we have chosen a rather small number of 
angular bins. The largest angular distance is determined by 
the survey geometry (Sect. 12. 1 ^ - for a patch geometry, we 
choose this to be slightly less than the patch diameter R. 
This is to avoid very large separations with small numbers 
of galaxy pairs, which caused our results to be instable when 
we varied the binning. 

The joint 2PCF data vector £tot has 38 entries, corre- 
sponding to the two times 19 angular bins for £+ and 
We calculate the aperture mass (Mf p ) for 14 different an- 
gular bins. The calculation of the covariance of (Mf p ) as 
a weighted sum of the 2PCF covariances can be biased for 
small aperture radii, when for a given binning of the 2PCF 
only a small number of terms contribute to the sum. We 
take into account a smaller number of aperture radii than 
separations for the 2PCF. 

3.2 Cosmological model and parameters 

We include m — 7 parameters in our analysis. These are f2 m , 
the power spectrum normalisation erg, the spectral index of 
the initial scalar fluctuations n a and V which determines the 
shape of the power spectrum. In addition, the cosmological 
constant A is included in our analysis. The redshift distribu- 
tion of source galaxies is assumed to be (Smail et al. 1995) 

p ( z )dz = ^— ( — ) 2 e- (z/zo) "dz. (17) 

The parameters zq and /3 are treated as free parameters to 
be determined from the data. 

Our reference model is a flat ACDM Universe with 
fi m = 0.3, r = 0.21, n s — 1 and erg = 1. The source redshift 
distribution is given by zo — 1 and (3 = 1.5, corresponding 
to a mean redshift of ~ 1.5. The CDM transfer function is 
taken from Bardeen et al. (1986), the non-linear power spec- 
trum is calculated using the fitting formulae of Peacock & 
Dodds (1996). 

3.3 KL eigenmode analysis 

3.3.1 Error bars 

For the case of a constant mean ('Sect. l2~4.1ll . we plot the 

MVB A Pa = (F QQ )- 1/2 = (2/ELi A fc) 1/2 in Fig. Has a 
function of the compression mode number h, for both single 
and joint parameter estimation. Typically, the asymptotic 
limit is reached at a lower mode number in case of £ + com- 
pared to On the other hand, the joint use of £ + and £~ 



reduces the error bars by a factor of nearly two in most cases. 
The saturation limit is approximately the same for each of 
the individual cosmological parameters. Irrespective of the 
survey strategy, approximately the first half of eigenmodes 
contain virtually all information about each individual pa- 
rameter, thus a compression factor of nearly two is possible 
without increasing the error. 

The error bars attained by both the individual and the 
joint analysis of the 2PCF £+ and £_ are typically tighter 
compared to the (Mf p )-constraints for a given survey geom- 
etry. 

We also conduct a joint parameter estimation for the 
case of a constant covariance (Sect. I2.4.2H . The complete 
information of all seven cosmological parameters is encoded 
in the m = 7 individual eigenmodes. However, the result 
of the SVD shows, that for f2 m ,o"8,.zo and (5 already the 
first singular mode carries basically all of the information 
on these parameters. For F, n s and A, the saturation of the 
error is reached after two modes. Apparently, the first mode 
picks up most of the information about the first group of 
highly-degenerate parameters, the second one completes the 
information about the parameters from the second group. 
This picture is consistent with the correlation matrix of the 
parameters, as discussed in Sect. 13.331 

3.3.2 Window functions for the compressed eigenmodes 

In Fig. |3 we plot the rank-ordered window functions Wi(£) 
associated with f2 m for the constant mean case. As an ex- 
ample, the first three window functions containing most of 
the information are compared with two higher order modes 
which contain less or negligible information. The survey 
strategy used in this analysis is a (30, 100')-patch geome- 
try (see Sect. EH- 

Because the window functions for £f are linear com- 
binations of Jo(£0j), there is always an extended tail which 
takes contributions from very large scales. The window func- 
tions for (Map) are very locallized and each eigenmode sam- 
ples only a small ^-region. For all statistics, small scales 
are noise dominated and contribute with only small ampli- 
tude to higher-order, low signal-to-noise KL modes, see also 
Fig. El 

The fact that the information content is larger for small 
£ and decreases with £ can be understood when consider- 
ing the simple case in which the covariance matrices and 
its derivatives are diagonal. Then, the eigenproblem matrix 
A = L'CqL - * H13H is also diagonal with Akk = (C :a )kk/Ckk 
being the eigenvalues. The diagonals of both C and C a are 
decreasing functions of the angular scale (except for small 
bumps due to geometrical survey effects). Since the decrease 
of C a is in general shallower than C, the largest eigenvalues 
occur on the largest angular scales. 

For the constant covariance case (Fig. the window 
functions Wo are much broader since there is only one eigen- 
mode containing all information and taking contributions 
from all angular scales. The Wo peak at a median angular 
scale where the signal dominates over the noise both from 
the intrinsic ellipticity dispersion of galaxies at small angu- 
lar scales and from the finite sky coverage at larger angular 
scales. 

The strong degeneracy among various parameters is re- 
flected in the behaviour of their associated window func- 
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Figure 3. Window functions Wi(£) multiplied with the corre- 
sponding eigenvalues A;, associated with the individual determi- 
nation of the parameter Q m , for the constant mean case. The first 
three eigenmodes, representing the three highest signal-to-noise 
modes, and for comparison, two exemplary higher order eigen- 
modes (i = 14, 19 for and £ t ot; i = 10, 14 for (Mf p )) arc 
plotted. The four panels correspond to the four different shear 
statistics as indicated. The specific survey strategy used in this 
case is a (30, 100') patch geometry (see Sect. |2~T1 . Note that the 
window functions for the two higher order modes have been dis- 
placed from zero to increase the visibility. 



tions: for all four statistics, degenerate parameters have very 
similar curves (see next section for a discussion of the differ- 
ent near-degeneracies). The F- and n s -filter functions have 
a zero-transition - the low ^-plateau has opposite sign than 
the peak reflecting the sensitivity of these parameters to a 
tilt in the power spectrum. 

Figure |H] displays for the dominant contibution to the 
Fisher matrix (second term in ea . 11011 how individual Fourier 
modes are sampled and combined in an "optimal" way to 
constrain cosmological parameters as indicated. It shows on 
which scales the convergence power spectrum has the largest 
influence on the determination of cosmological parameters, 
depending on the survey strategy (see also Sect. l.'-i. 3~4*)l . In a 
forthcoming paper, analogous analyses based on the present 
work will be performed for third-order statistics of cosmic 
shear, allowing one to quantify the sampling of the con- 
vergence bispectrum to extract cosmological information. 
Moreover, the scales of interest will be studied for a com- 
bined analysis of the power and the bispectrum which will 
reduce the near-degeneracies between parameters. 



3.3.3 Correlation between cosmological parameters 

In Fig. we plot the correlation coefficient of the Fisher 
matrix r a fj = F a p / ' (FaaFpp) 1 ^ 2 , which quantifies the level 
of degeneracy between two parameters, as a function of the 



Figure 4. The dominant first window functions \\W\(t) for the 
constant mean case, associated with various parameters as indi- 
cated in the figure. The four panels correspond to the four differ- 
ent shear statistics as indicated. The survey strategy is (30, 100') 
(see Sect. I2TI . For visual clarity, some of the curves are shifted 
vertically. 
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Figure 5. Window functions cnWi for the joint analysis of all 
parameters using SVD in the case of constant mean. The first four 
singular modes multiplied by the corresponding singular values 
are displayed. 
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Figure 6. The window function Wo for the constant covariance 
case associated with various parameters as indicated in the figure. 
The four panels correspond to the four different shear statistics 
as indicated. 
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Figure 7. The correlation of H m (upper panels) and Y (lower 
panels) with the other parameters, as a function of the mode 
number n' for the joint parameter estimation involving all pa- 
rameters. The left and the right panels correspond to £tot and 
(M 2 p ), respectively. 





Figure 8. Ap a as function of the eigenmode number ft, for the 
survey geometries (30,100') (solid lines) and 75 ■ 26' 2 (dashed 
lines). The single parameter case of Q m ,crg and Y is displayed. 



compression mode number h. For very small mode num- 
bers, the Fisher matrix becomes singular and all parameters 
are completely degenerate. Clearly, there exist two groups 
of highly-degenerate parameters (f2 m , erg, zo, P, Oa.) and (r, 
n s ). For parameters from different groups, the correlation 
decreases, even for mode numbers h where the diagonal ele- 
ment has reached the plateau. Thus, redundant modes which 
do not carry any information regarding the Fisher diagonal 
are nevertheless important and help to reduce parameter 
degeneracies. 

We also calculate the correlation coefficient r a p for the 
constant covariance case, where the SVD gives us only m — 7 
singular modes. The correlation coefficient forms a plateau 
when as few as two of the singular modes are included, as 
it is the case for the diagonal elements of the Fisher matrix 
(Sect. 1X3. H . However, even though F does not change much 
when three modes and more are added, this is not true for 
the inverse Fisher matrix since F is quite ill-conditioned and 
small changes have a large impact on F _1 . For a numerically 
stable determination the MVBs from the inverse Fisher ma- 
trix, one needs as many modes as number of parameters. 



3.3.4 Survey strategy 

We plot Ap a as a function of the eigenmode number in 
Fig. |H| In contrast fo Fig. [5] where we compared the single 
and the joint (SVD) parameter case for one survey strategy, 
we emphasis here on the single parameter estimation and 
compare the resulting errors for different survey geometries. 
All patch geometry strategies as well as the survey consist- 
ing of 12 uncorrelated 65' 2 -fields yield very similar results, 
the case (30, 100') is shown as a representative. On the con- 
trary, the MVBs from the 75-26' 2 -geometry are significantly 
higher with few exceptions. This latter survey does not sam- 
ple medium and large scales in contrast to the patches and 
the 12 ■ 65' 2 case. The small cosmic variance corresponding 
to the 75 independent lines of sight can not compensate for 
the lack of large-scale information. 

The effect of survey strategy on the window functions is 
a shift in £, depending on the scales that are sampled by the 
survey (compare Fig.lBlwith llOl and Fig.lolwith lllll . For the 
constant mean case, the uncorrelated images-surveys shift 
the peaks of the (M 2 p )-window functions corresponding to 
f2 m and the other parameters from this degeneracy group 
(Sect. 13.3.31 towards higher I whereas for F and n s , the 



10 Kilbinger & Munshi 



( M od) 




5 


10 


SVD 


(M 2 p ) 


i < i i i 
/\ 

. i V , \ 


i i /I 




5 


10 





1 I , s 

-1 ' 1 - 


A A, A, 




-■' 1 / I'M 


/ A t , 1 I" 


I 


A yt 



5 10 

n 

(30,100') 

(30,140') 

75.26' 

12.65' 



Figure 9. The peak position of the n th (Mf p )-windows as a func- 
tion of the corresponding eigenmode number n, for single analy- 
sis of the parameters O m (upper left panel) and T (upper right), 
and for the joint estimation (lower panel). £* P eak = 5/^pcak = 
17'(1000/^ pca ^) is the corresponding peak in real space. 



peaks seem to be randomly distributed, see Fig. [5] For the 
latter two parameters, there is no preference of large scales, 
apparently all scales contribute with the same importance 
to the Fisher matrix. 

Figure TPH shows the window functions Wo for the joint 
parameter analysis, corresponding to the dominant constant 
covariance case. From this figure, it is clear that only the 
first two modes carry significant information about cosmol- 
ogy. The shapes of the two functions is complementary, for 
example, the second one has a zero transition at roughly the 
peak position of the first one. The reason is that there are 
basically two groups of parameters entering the convergence 
power spectrum, and the information content is described 
by two independent singular modes. 

For various survey strategies, we calculate the simul- 
taneous MVBs for f2 m ,<T8,r and n s using ©. The patch 
geometries yield similar errors which are smaller than those 
from the uncorrelated images-surveys, see Table 13.3.41 In 
particular for (Mf p ), the lack of large scales of the 75 • 26' 2 
case results in a large uncertainty for fi m and ag. This con- 
firms a result found by Jain & Seljak (1997), stating that in 
order to break the f2 m -<78-degeneracy, the shear correlation 
on large (linear) scales has to be added to the measurement 
of small (non-linear) scale correlation. 



3.3.5 Noise levels 

We determine the dependence of the Fisher matrix for differ- 
ent survey characteristics and noise levels. The calculations 
are done for a (30, 100')-patch survey geometry. We vary the 
survey area A between 1.4 and 14 square degrees, by adding 
more and more independent, uncorrelated patches, each con- 
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Figure 10. The first four window functions aiWi for the joint 
parameter and constant mean case, corresponding to a 75 • 26' 
survey. See Fig. [F] for an explanation of the curves, and for com- 
parison with a (30, 100') geometry. 
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Figure 11. The constant covariance window functions Wo for 
the seven cosmological parameters, corresponding to a 75 ■ 26' 
survey. See Fig. |H] for an explanation of the different curves, and 
for comparison with a (30, 100') geometry. 



Table 1. Simultaneous MVBs for various parameters and survey 
strategies. 
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Figure 12. The first four window functions multiplied by the 
singular values, CjWi,i = 1 ... 4, for the joint analysis of constant 
covariance, corresponding to two survey strategies as indicated. 



taining 30 images. For the variation of the other parame- 
ters (source ellipticity dispersion o> , number of background 
galaxies n ga i, number of bins n and redshift parameter zo), 
a single patch was used to calculate the Fisher matrix. In 
the constant covariance case, we write the MVB as 



Ap a = F c 



'1/2 



= k 



A 



10 sq deg 

/ "gal 



0.3 



V20J 



■■ , I -),» / ■" (18) 
, 30 arcmm / V 20 J 

where the constant k and the power-law indices fj,, v, n and 
A are given in Table [5] for each parameter a. 

Note that although the MVB for (Aff p > is smaller than 
the one for the 2PCF, this is not true when more than one 
parameter is considered. In that case, the MVB is given by 
the inverse of (a submatrix of) the Fisher matrix. The off- 
diagonal terms are typically larger for (Mf p ) which decreases 
the diagonal of the inverse matrix. 

If the mean is constant, the Fisher matrix is indepen- 
dent of the area of the survey. Except for the binning (see 
next section), the sensitivity of the other survey character- 
istics is weaker than for the constant covariance case. For 
our fiducial values of a e = 0.3, n ga i = 30arcmin~ 2 ,n = 20 
and zo = 1, the survey area where both terms in il l (It are 
equal varies between 0.06 (for n s ) and 0.4 (H m ) square de- 
grees in the case of £tot- For the aperture mass statistics, 
this area of equal contribution of both terms to the Fisher 
matrix is roughly a factor of 4 smaller. Thus, only for very 
small survey areas the constant mean term is important. 

Noise at small angular scales is dominated by Poisson 
noise coming from the intrinsic ellipticity dispersion of the 
galaxies <r e , and the finite number density of galaxies n ga i. 
The effect of both error sources is shown in Fig. [T31 The 



Table 2. The coefficients describing the MVB for the constant 
covariance case as function of various survey characteristics, for 
different cosmological parameters, see 1181 . 
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0.0074 


-0.77 


0.40 


-0.03 


0.52 


r 
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information content of each eigenmode decreases with de- 
creasing n ga i and increasing a e . The asymptotic plateau is 
reached for marginal smaller mode numbers in the case of 
higher noise level. The resulting decrease in accuracy with 
which a given parameter is determined can be inferred from 
the asymptotic values. 



3. 3. 6 Binning 

We vary the number of angular bins n of the two-point cor- 
relation functions between 20 and 50, in order to quantify 
the effect on our results. For the constant mean case, the 
power-law decrease of the Fisher matrix diagonal are unaf- 
fected by the binning. We found the height of the plateau 
to weakly depend on the bin number, as n~ v where r\ is 
different for different statistics and parameters, and ranges 
between 0.18 and 0.3. 

Note that the compression factor is somehow arbitrary 
since the number of bins can be chosen freely. However, 
within the (realistic) range that we tested, it is close to two. 
If wc further increase the number of bins, the compression 
will of course increase; no new information will be added 
since the correlation between adjacent bins will approach 
unity for very small separations. The weak dependence of 
the plateau on n shows that we have not yet reached that 
saturation regime. If one considers the observed galaxy el- 
lipticities as input data vector, the compression factor is of 
course much larger. 

The Fisher matrix in case of constant covariance is only 
very weakly dependent on the bin number, as shown in the 
previous section (see Tabled- The diagonal elements of the 
Fisher matrix for both cases (constant mean/covariance) are 
shown in Fig. 1141 In both the cases, for constant covariance 
and constant mean, the amplitude of the window functions 
increases with the number of bins. The shape is not affected 
by the binning. 
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Figure 13. Ap a = (Fan) -1 ' 2 as a function of the number 
of modes n as in 1141 . for different noise levels. Upper panels: 
n gal = 25arcmin — 2 (solid lines), n ga i = 50arcmin -2 (dashed 
lines). Lower panels: a e = 0.2 (solid lines), cr e = 0.4 (dashed 
lines). The left panels corresponds to the joint correlation func- 
tion £tot whereas in the right panels results for (Af 2 p ) are shown. 
The plots are done for the parameters f2 m , erg , T and n s . The sur- 
vey geometry is a single (30, 100') patch. 



4 DISCUSSION 

We study the Karhunen-Loeve generalised eigenmode prob- 
lem in the context of weak lensing surveys. Four combina- 
tions of second-order shear statistics estimated at various 
angular scales are used as input data vectors - these are the 
2PCFs £±, their combination £ t ot and the aperture mass 
dispersion (Mf p ). Different realistic survey geometries are 
considered to study the impact of noise and finite sky cov- 
erage for these estimators. We examine the dependence of 
both the mean and the covariance of these estimators on var- 
ious cosmological parameters and the source galaxy redshift 
distributions. A Fisher matrix analysis is presented to de- 
termine errors associated with the estimation of various cos- 
mological parameters and their cross-correlation with survey 
parameters. 

We consider two different scenarios. In one case, the 
mean of second-order shear estimators is used to constrain 
the cosmological parameters while the covariance is con- 
stant and independent of cosmology. The second case uses 
the covariance of these estimators to constrain cosmolog- 
ical parameters assuming a constant mean. We study the 
information content of KL eigenmodes in both cases for var- 
ious parameters. For the constant mean case, there are sev- 
eral eigenmodes among which the information is distributed; 
however, there is only one eigenmode associated with the 
constant covariance case which contains the complete in- 
formation. The resulting error bars are anti-proportional to 
the square root of the survey area for the constant covari- 
ance case, and independent of the survey area if the mean 
is constant. Thus, for reasonable sky coverage (more than 



Figure 14. Effect of the bin number n on the Fisher matrix in 
the constant mean case for various estimators as indicated in the 
panels. 

about 0.4 square degrees) the first case dominates over the 
second and the bulk part of the cosmological information is 
collected in only one mode. 

From the results of our KL eigenmode analysis we find 
that, starting from n binned values of second-order shear 
statistics, a compression factor of almost two can be achieved 
in most cases without information loss, if the cosmological 
parameters are determined using the covariance. KL analysis 
provides rank-ordered, uncorrelated eigenmodes and corre- 
sponding window functions in the case of single parameter 
estimation. These sets of window functions provide clean- 
est measures of the projected power spectrum for a spe- 
cific survey strategy. The first eigenmodes contain most of 
the information that one can extract and directly use for 
maximum-likelihood studies to constrain cosmological pa- 
rameters. Typically, the error bars on individual parameters 
using the covariance are smallest when the joint correlation 
function £ to t is used. gives better constraints than £_ 
which is in turn better than the aperture mass statistics 
(Mf p >. 

We provide both an independent study of different pa- 
rameters as well as a joint analysis of all parameters for each 
of these two cases (constant mean and constant covariance). 
The parameters which are near-degenerate have similar KL 
eigenfunctions. Among the five cosmological parameters we 
have considered, (f2 m , as) and (r, n s ) show similar levels 
of degeneracy and have very similar eigenmodes. Moreover, 
the parameters characterising the source galaxy redshift dis- 
tribution (zo and 13, see 11711 ) are degenerate to a high level 
with Q, m and a$. 

It may be useful to note that additional constraints 
may be put by including higher-order statistics which will 
be left for further studies. Moreover, since the number of 
observables increases with increasing order (i.e. the shear 
three-point correlation function has eight independent com- 
ponents and depend on triangle configurations, Schneider & 
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Lombardi 2003), data compression can be useful and even 
necessary. 

The eigenmodes associated with fi m and er 8 mainly fo- 
cus on larger angular scales where cosmological information 
is not contaminated by noise on small scales. On the other 
hand, F and n s measure the shape of the projected matter 
power spectrum which is easier to determine when also in- 
formation on smaller angular scales is available. Therefore 
the eigenmode windows associated with these two parame- 
ters tend to take more contributions from smaller scales. In 
contrast, the eigenmodes associated with the constant co- 
variance case take contributions from virtually all angular 
scales probed by the survey and reaches a maxima roughly 
at medium angular scales where the signal dominates over 
both shot noise and scatter due to finite sky coverage. 

Higher order KL modes do not carry any useful infor- 
mation and approximately half of the modes contain all of 
the information content of the Fisher matrix for a constant 
mean. However, the off-diagonal terms of the Fisher matrix 
which encode the cross-correlation among various parame- 
ters are also a function of the number of modes we include 
to reconstruct the Fisher matrix. We found that these terms 
keep evolving even after the diagonal terms have already 
reached a saturation limit. For very small numbers of modes 
the reconstructed Fisher matrix becomes singular, since the 
information sampled in only the few eigenmodes is too lit- 
tle in order to put constraints on more than one parameter, 
thus different cosmological parameters become completely 
degenerate. 

The joint analysis of all seven parameters being esti- 
mated from the mean only (constant covariance) shows that 
only two eigenmodes (out of seven) are needed to constrain 
the parameters and to lift their near-degeneracies as far as 
possible. 

In our analysis involving KL eigenmodes we have not 
included and modelled any systematic measurement errors. 
We have assumed that the errors are dominated by Pois- 
son noise and intrinsic ellipticity noise at small scales and 
cosmic variance at larger angular scales. The covariance 
of the second-order shear statistics include all these noise 
sources, which are exact if the shear is a Gaussian field. 
Non-Gaussianity leads to an under-estimation of the noise 
on angular scales between 1 and 10 arc minutes. 

Although the exact implementation of the survey ge- 
ometry makes our approach more realistic, we still ignore 
B-mode contamination due to systematic residuals. This is 
mainly due to the absence of a specific model for the B-mode 
power spectrum in the case of systematic measurement er- 
rors. Such issues will have to be included in future studies. 
Although effects like intrinsic galaxy alignment and source 
clustering, which also cause a B-mode, can be and have been 
modelled (Crittenden et al. 2001; Croft & Metzler 2001; Hey- 
mans et al. 2004; Schneider et al. 2002), their contribution 
to the shear signal is secondary. 

Incorporation of (photometric) redshift information will 
in general improve the accuracy with which cosmological pa- 
rameters can be recovered. However, we have not divided 
the source galaxies into redshift bins, because of the corre- 
sponding complication in the calculation of the covariance. 
However, this is possible in principle, e.g. using correlated 
Gaussian fields, see e.g. Simon, King & Schneider 2004. The 
inclusion of non- Gaussian terms in the Fisher matrix anal- 



ysis (Taylor & Watts 2001) will invariably help us to break 
some of the parameter degeneracies which we have stud- 
ied here. However, the construction of covariance matrices 
associated with various third-order estimators tend to be 
very cumbersome. Although certain analytical calculations 
were made recently under simple assumptions (Munshi & 
Coles 2003; Munshi & Valageas 2005), more detailed analy- 
ses which incorporate realistic survey geometries remain to 
be done. 

We incorporated the non-linear model of the power 
spectrum according to Peacock & Dodds 1996. In a forth- 
coming work, an improved version based on the halo model, 
e.g. halofit (Smith et al. 2003) will be used. In addition, 
baryonic physics, non-zero neutrino masses and a varying 
dark energy equation-of-state parameter will be taken into 
account. 

A simple-minded signal-to-noise eigenmode analysis (an 
interesting special case of KL eigenmode analysis, see TTH) 
is difficult to perform for the case of the covariance matri- 
ces we have considered. This is because of the difficulty in 
separating the signal and noise contributions in the covari- 
ance matrices. Where as the V terms (cosmic variance) are 
pure signal and the D terms (shot noise) is pure noise, it is 
difficult to separate the mixed terms M with such clarity. 
Increasing the noise contribution either by increasing the 
intrinsic ellipticity dispersion of galaxies o> or by decreas- 
ing the number density of source galaxies n ga i results in a 
decrease in the information content of each of the KL eigen- 
modes. In addition the number of eigenmodes which contain 
useful information reduces, too. 
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